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ABSTRACT 


The buckling aotion of hollow cylinders, subjected to internal 
pressure, was investigated under various conditions of axial load in the 
oylinder walls. The theory of Euler was found to apply where the modulus 
of elasticity and second moment of cross-sectional area are determined 
solely by the material and cross-section of the oylinder, and the load 


is the total axial load on the column, 
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SYMBOLS 


Cross-seotional area of the oylinder (in2). 


Cross-sectional area of the piston (ino) 


Lateral deflection of the mid-point of the column 
[άν τη. κ 


Modulus of elasticity of the material of the 
cylinder (psi). 


Second moment of cross-sectional area about a 
centroidal axis (in.%). 


Length of the column (in.). 
Bending moment (in.-lbs.). 

Axial Compressive load (lbs.). 
Critical load (lbs.). 

Internal pressure (psi.). 

Critical internal pressure (psi.). 
Radius of ourvature (in.). 


Radius of gyration (ina). 





I Introduction 
1. Introduction. 
A long cylindrical shape subjected to a compressive losd will fail 
rtv sudden ben lin, provided the membor is slender. This is the type of 
loading investi ated my [ular [1] who doternined a "critical load" such 
that an load treater than this will cause the column to bend. In his 
anal sis culer considered the column to be solid and subjected to an axial 


campressive load applied at the ends. 





| Urs 1. Sohomatis tepreserntation of Column 


assuñte lindar fillad with а fühid capa*le o7 be?n? 


o l uq r (resskre (-i7. 1). m3 ond of the cylinder is ricidlv sealed 


L) 


end the os#er “5 sealed by means of a Srictionlsss viston. Let tne fluid 
ua vlacei infer pressure: by means o" an axial compressiva load P applied 
at $no ende. Synbelicallvy the 5heory o? wuler states: і? p is the 

critical load; 5 the modulus of elasticit, o® tie material; I the second 


moment of the cross-sectional area about a cer.roidal axis; and 4 the 





the lenrth of the column in the hinged-end oconfiruration; then 
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A free body diagram of the column of Fig. 1 is identical to the free body 





diagram of any column with the same external restraints. Applying Euler's 
equation to the colunn of Fig. l results in: The load being the product 
of the piston aren A. p and the internal pressure 4. Ihe modulus of 
elasticity and the second moment of the oross-sectional arsa are deter- 
mined solely by the material and the cross-section of the cylinder. 
Writing Euler's equation in terms of a hollow cylinder filled with fluid 


results in 





| „ untl 
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In the undefleoted condition the entire axial load of the column pf 
Figo 1 is carried by the fluid. Thus we are led to the interesting 
paradox of a column, which in the undefleoted condition sustains no axial 
load, buckling when the internal pressure reaches a certain oritiosl value. 

Interest was fooused on this problem when a long cylinder used in the 
hydraulic system of a submarine was observed to buckle when the piston 
approached the end of its stroke. With the increase in operating pressure 
of hydraulio plants and the associated decrease in size of component 
parts, it 1s oonoeived that in some applications the buokling action of 
oylinders subjected to internal pressures will be important. 
2„ Object and Scope 

The object of this investigation was to verify the applicability of 
the theory of Euler to hollow cylinders subjected to internal pressure. 


It was previously shown that 





πα: 
Ф.) = A, γ᾽ 


This equation states that, for a given length, the internal pressure which 
will cause columns of the same material to buckle is proportional to the 
quotient of the second moment of cross-sectional area by the area of the 
piston. 

The axial load carried by the cylinder is determined by the relation 
between the area of the piston and the area of the inside of the cylinder. 
The area of the inside of the cylinder does not appear in the equation. 
Therefore it may be stated that, theoretically, the critical pressure is 
independent of the "load-ratio." (Load-ratio is defined as the ratio of 
the axial load carried by the cylinder to the total axial load on the 
column. ) 

To verify this conolusion it was planned to test three sets of 
columns. In the first set the inner area of the cylinder was less than 
that of the piston (a positive load-ratio). In the second set the areas 
were equal (a zero load-ratio). In the third set the inner area of the 
cylinder was greater than that of the piston (a negative load-ratio). 

To check the validity of the theory through a range of pressures, 
the lengths of the column were changed so that slenderness-ratios from 
100 to 200 were obtained. The effeot of these variables was studied 
experimentally and the results correlated with theory. 
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II Theory 

Consider a straight cylindrical tube of homogeneous material, wnich 
obeys Hooke's Law, filled with a fluid capable of being placed under 
pressure, One end of the tube is rigidly sealed and the other end is 
sealed by means of a frictionless piston. The two ends are subjected to 
an externally applied axial load. The ends are restrained against lateral 
movement but have rotational freedom. The free body diagram of the column 
(Fige 1) is identical in form to the one considered in Euler's classio 
worke 

From Euler's solution it is known that under an axial compressive 
loading of magnitude P, applied at the oentroids of the end cross-sections, 
the tube will remain straight until P reaches a certain critical value. 
As P exceeds this oritical value the tube will begin to bend laterally in 
a plane containing the axis of the tube, If the tube has a circular cross- 
section this plane will not be known in advance. The lateral bending 
increases rapidly with P, so that it is practical to consider that the 
critical value of P represents the maximum load that the column can 
support. In order to evaluate this critical load we choose a set of 
cartesian coordinates, with the x-axis the locus of centroids of the oross- 
sections of the undeformed tube, the xy-plane containing the axis of the 
deflected tube, and the z-axis perpendicular to the x and y axes. The 
origin is at one end of the tube. If we consider the column in a deflected 
position (Fig. 2) it is olear that the magnitude of the bending moment at 
a Pert ioular cross-section is Py where the y-ooordinate is that of the 


centroid of the section. 
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The condition of particular interest fa the ara in which the column first 


starts to deflect. At this time the 


Slope is snall and nerlectin~ it 
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will produce a neglizible error and creatly sinpli^y the solution. 
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Thus we have 2 


|. a d y 
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o Ón—— А 


Zn dx” 


Employing a consistent set of algebraio signs we are led to the differen- 


tial equation D 
Ae 
Ti У = (О. 
ах El 
The general solution of this equation is 
Y = А с RX +B am Rx, where 
R 
REI 
The end oonditions are satisfied if 
and either: B =O 


EE on RA =0 
If D:O then ZO and this is the solution prior to buokling. 


If _0mRkl:O then 
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mnt El 
or P Ξ JT where m is an integer. 


The solution of principal interest is the one obtained for m» | , 


from which 2 
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cr Ё? 

(The solution obtained for uO is that of the undeformed column. 


Solutions obtained for values of m greater than one require that the 
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Fizure 3. Schematic Representation of the “xperimental Set-up 


ZA 


Fiure 3 is a schematic representation o: the experimental set-up. 
ir the piston is frictionless then P may be replaced ty its equivalent 


in terms of interral pressure RA, . The fluid can not sustain shear 


so that the modulus of elasticity; and the second moment of oross-sectional 





area depend only on the material and cross-section of the cylinder. 


Placing the equation in terms of the internal pressure we obtain 


(4). - BEA | 


Where A. is the area of the cross-section, and © is the radius of 


gyration of the oross-section of the cylinder, 





III Description of Apparatus 

1. Selection of Material 

The bucklin5 load, and thus the internal pressure producing buckling, 
is directly proportional to the modulus of elasticity. Aluminum alloy 
has the lowest modulus of elasticity amonz the readily available tubing 
materials and was therefore chosen. 
2e Choice of Geometrical Properties 

The principal limitation placed on the geometric properties was that 
the tube buckle before the hoop tension reached the yield point. A 
slenderness-ratio of 100, which fulfills this requirement, was selected 
as the minimum to be tested. | 

Establishing a minimum slenderness-ratio and a maximum internal 
pressure (see Section 5) dictates a maximum diameter-to-thickness ratio. 
Among the geometrically similar tubes which could be used outer diameters 
of 3/4, 1 and 1-1/4 inch were selected because they provide lengths which 
could be conveniently handled. A thickness of „035 inch was required for 
the 1-1/4 inch tube. This thickness also provided satisfactory charac- 
teristics for the other two sizes. 
š, End Fittings 

Figo 4 is a schematio diagram of the end fittings. The tube and 
end fittings are joined by extra fine threads, 56 threads per inch. The 
upper end fitting has a cylinder joined to it by threads. (Henceforth 
this will be considered one piece and called "the upper fitting.”) The 
upper portion of the piston, which contacts the horizontal surface of 
the testing machine, is a spherical surface (radius 2.75 inches). The 
eater of this surface is on the axis of the piston in the plane of the 
face of the upper fitting. (Face of a fitting is defined as that flat 


10 
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sarface tirourh w ¿ch the tube is t'roaded.) 





8) Lower b) Upper 
Fire A. "nd Fittings 


The lower fittin is similarly formed with the center of the 
spherical surface (radius 1.75 inches) on its axis in the plane of its 
face. The lower fitting has a passage, fitted with a 1/4 inch pipe 
nipple, to allow entry of the hydraulic fluid. 

Ihe effect of the spherical surfaces is to cause the column to act 
as if it were in tho hinged-end configuration with an effective lenrth 
equal to the distance between the faces, 

The piston and upp or fittings (Fir. 5-a) retain the hydraulio 
fluid and dictate the axial load transmitted to the tube. The piston 
is fitted with an O-rins seal to limit leakage. The same piston is used 
in testing the three different sizes of columns. This results in a 


positive load-ratio, a zero load-ratio or a nezative load-ratio; 
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b) Upper fitting 
with 3/4 inch 
tube 


a) Piston, and upper 
fitting with 1 | 
inch tube 


c) Upper fitting 
with 1-1/4 inch 


tube 








depending on whether the inner diameter of the tube is lesa than 
(Fig. 5-b), equal to (Fig. б-а), ог greater than (Fig. 5-0) the diameter 
of the piston. (The diameter of the piston is taken to be the inner 
diameter of the oylinder.) 
4, Deflection Measurement 

Lateral deflections were measured by means of two dial gages set at 
right angles to each other and mounted, at the mid-point of the column, 
in a plane normal to the axis of the undeformed columne 
5. Pressure Application 

The pressure was applied by a dead-weight gage testing machine, 
using hydraulic fluid. The apparatus was designed to oarry pressures to 
2000 psi. It was felt that pressures greater than this would complicate 
the sealing of the threaded section and require a more complex O-ring 
designe 

The fluid is carried from the gage tester to the lower fitting 
through 10 feet of hose having an inner diameter of 3/16 inch and an 
outer diameter of 33/64 inch. The hose consists of a synthetic rubber 
inner tube, a single wire braid reinforcement, and a synthetic rubber 
cover. 

Fig. 6 shows & photograph of the apparatus assembled and ready for 


test. 
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Figure 6. Photograph of the Apparatus 





IV Experimental Procedure 
In derivin; the theory it was assumed that the piston was friction- 
less. The experimental procedure was devised to eliminate the effects 
of friotion between the piston and oylinder. The first method attempted 
was to apply the pressure, then turn the piston within the cylinder. 
This method appeared to eliminate the friotion, but the resulting lateral 
movement of the column made the lateral deflection readings unreliable. 
The method finally used to eliminate the friotional effect was to apply 
the pressure and then adjust the testing machine to a load equal to the 
product of the piston area by the applied pressure. The column was 
aligned in a vertical position by using the plumb bubbles on a 28-inch 
carpenter's level. 
The sequence of operation was as follows: 
1) The column was placed in a vertical position by 
means of the level; 
2) a fraction of the oritioal pressure was applied 
by loading the gage tester; 
3) the load on the testing machine was adjusted to 
that load which it would oarry if the piston were 
friotionless; 
4) the lateral deflections were read from the dial gages; 
5) the piston displacement was checked to insure that it 
had not reached the end of its travel; 
6) the pressure applied by the gage tester was inoreased, 
Steps 3) through 6) were repeated until the defleotion read in Step 4) 


reached the limit of the dial gage. The oritioal pressure was 
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ucrested by Southwell (Appendix I, Sect ion 1). 
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V Results and Conolusions 
1. Tabular Results : 

The theoretical and experimental values of critical pressure and 
unit load (unit load is the quotient of the total axial load on the 
column by the cross-sectional area of the oylinder) are summarized in 
Table 1. Additional data appears in Appendix I, Section 4. 

2. Graphical Results 

The theoretical curve of unit load versus slenderness-ratio, with 

experimental points plotted, is given in Fig. 7. 


3. Conclusions 


The experimental data indicate that the equation of Euler 


D - ΠΕΙ 
Z 


is applicable to hollow cylindrical shapes subjected to internal 





pressure, where: P is the total axial load on the column; and the 
modulus of elasticity and second moment of cross-sectional area depend 


solely on the material and oross-section of the cylinder. 
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Table 1 


Column Slenderness Critical Pressure (psi) Unit Load % Dif- 
Number ratio Theoretical Experimental Theoretical Experimental ference 
3/4 inch tube 
1 105.1 1092 1074 9460 9500 -1,7 
2 105.1 1093 1091 9460 9450 -0.2 
3 120.2 835 820 7220 7090 - 
= 120.2 835 819 7220 7080 -1.9 
5 147.8 999 545 4790 4690 -2.2 
6 147.8 553 543 4790 4700 -1.8 
7 202.6 296 296 2550 2550 ο 
8 202.6 297 292 2550 2510 -1,7 
l inoh tube 
9 100.5 1643 1625 10360 10250 -l.1 
10 100.5 1646 1649 10360 10380 20.2 
ll 118.7 1172 1189 7430 7540 21.4 
12 115.1 1190 1194 7500 7520 5.0.3 
13 142.2 816 800 5180 5070 -2.0 
14 142.5 816 811 5170 5140 -0.6 
15 185.6 474 465 5040 2980 -1.9 
16 185.5 475 474 5040 5030 =005 
1-1/4 inoh tube 
17 101,5 2068 2012 10160 9890 - 2 7 
18 101.5 2069 2024 10160 9940 - 2. 2 
20 118,0 1508 1503 7510 7480 -0.35 
21 136.5 1138 1120 5650 5540 -1.6 
22 136.5 1136 1118 5630 5950 -1.6 
23 152.4 917 Un 4510 4510 0 
24 152.4 219 DD 4510 4510 -0.2 
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APPENDIX I 
ANALYSIS OF EXPERIMENTAL DATA 

1. Southwell's Method 

In the ideal case of column buckling there would be no lateral 
deflection up to the critioal value of the load. At the oritical load 
the column would then begin to deflect. Owing to various kinds of imper- 
feotions, such as some unavoidable initial ourvature of the colum, 
eccentricity in application of the load and non-homogeneity of the 
material, the column begins to deflect with the beginning of loading and 
usually fails before Euler's load is reached. 

A useful method of determining the critical load from the test data 
was proposed by Re V. Southwell [ο]. Assume that the deflection of a 
column, under a load that is below the critical value, is due to initial 
curvature and eccentricity of load application. If the deflection is 
expressed as a trigonometric series, when the load approaches the critical 
value the first term in the series predominates. 

If the deflection d is represented by the first term of the series, 


the following linear relation exists 


gle [^ š O 


where Q is a constant. A cartesian plot of d versus £ yields a curve 
whose slope asymptotically approaches the critical load. 
The experimental data of Column 11 are plotted in Figs. 8 through 


12, The nature of tho curve and its asymptotic behavior can be observed, 
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2o Sample Caloulations 

The critical load used in the theoretical derivation is the product 
of the internal pressure and the piston area. In the experimental 
procedure it was advantageous to consider the internal pressure to be 


the principal variable. For this reason the data are reduced in terms of 


the internal vressure. 
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Column 11 Run 1 


Gage 1 Gare 2 

+ d d/p, d a/a 
200 209 5 «025 
400 3.5 15 .052 
600 45 22 0037 
700 5 52 .046 
800 6 42 .054 
850 6 49 .058 
900 6.5 65 .070 
950 6.0 T2 «076 
975 6 84 «086 
1000 6 100 9100 
1025 6 108 e105 
1050 + 155 ‚127 
1075 4,9 160 «149 
1100 0,9 187 ¿170 
1125 -8 250 „222 


Gage l shows a possible variation of the data. At moderate 
pressures same term of the series other than the first predominates. 
As the critical pressure is approached the principal mode of 
distortion predominates. The data are not plotted since no trend 
can be determinede 


Data from Gage 2 are plotted in Fig. 8. 
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£ 
200 
400 
600 
700 
800 
850 
900 
950 
975 
1000 
1025 
1050 
1075 
1100 
1125 


1150 


Column 11 


Gage 1 
d 

1 
5.5 
705 
9,5 
15 
16.5 
18.5 
202.0 
2305 
265 
29.5 
52,5 
56,5 
59,5 
41.5 


47.5 


Data fran Gage 1 are 


Data from Gage 2 are 
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Run 


d/p. 
«0050 


o 0088 
. 0125 
«0156 
«0187 
©0194 
.0210 
.0257 
.0241 
©0260 
‚0288 
©0310 
|. 0540 
.0560 
00569 


00413 


Gage 2 


18 
24 
34 
41 
49 
64 
69 
83 
95 
110 
129 
155 
175 


250 


plotted in Fig. 9. 


plotted in Fig. 10. 


dh: 
.015 
.020 
.030 


| 054 


. 048 
«056 
«067 
.071 
.083 
.093 
9105 
«120 
„141 
«154 
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Column 11 Run 3 


Gage 1 Gage 2 

Ф d а/д d 3/4 
200 0 9 9020 
400 1 12.5 0031 
600 1 23.5 .059 
700 0 3200 «046 
800 0 43.5 «054 
850 0 9300 .063 
900 О 63.5 071 
950 -1 |.0015 73.5 . 078 
975 -3 „0030 88.5 . 091 
1000 -5 «0030 93.5 «094 
1025 -5 .0049 106.5 . 104 
1050 -8 .0076 138.5 0132 
1075 -11 ©0102 164.5 0153 
1100 -15.5 .0141 171.5 «156 
1125 -26 002351 236.5 „210 


Data from Gage 1 are plotted in Fig. ll. 
Here is an exemple of the first term of the serier 
predominating as the critical pressure is approached. 
In this case a trend is established. 


Data from Gage 2 are plotted in Fig. 12. 
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5. Criteria Used in Evaluating Critical Pressures. 
Let the points ottained fran the data be numbered "1" throu-h "n" 
in order of inoreasinz load. For a "good" data set the following 
generalization oan be made: 
1) The slope of the line joining points "n-2" and "n" will 
be approximately equal to the slope of the line joining 
points "n - 1" and "n," 

2) The best estimate of the asymptotic slope, without 
extrapolation, is the slope of the line joining points 
"n - 1" end "n," 

Table 2 lists the slopes obtained from the last three points of 
Figs. 8 through 12 

Table 2 


Slopes Obtained from Data Sets of Column 11 


Run 1 
Gage 1 Gage 2 
No curve 1251 
1210 
Run 2 
Gage 1 Gage 2 
1510 1235 
1362 1202 
Run 3 
Gage 1 Gage 2 
1163 1262 
1165 1202 
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The following criteria were set up to evaluate the information 
obtained fran the plots. Data sets in which the deflection ohanges 
direotion within the last three readings were discarded without plotting. 
Data sets in which the two slopes computed fram the last three points 
differed by more than five per cent were discarded. Data sets in whioh 
the two slopes agreed within five per cent were assigned a weight of 
unity. Data sets in which the two slopes were within one per cent were 
assiened a weight of 2. The weishted average of the slopes of the lines 
joining the last two points was taken as the critical pressure of the 
colimn. 

The above system of evaluation, though it does not prevent disoard- 
ing good information, makes it unlikel: that poor information will be 
included. 

4. Tabulation of geometrical properties and critical pressures. 

The geometric properties, oritical pressures, aid average critical 
pressure for each column are given in Table 3. Method of determining 
theoretical values used in Table 3 may be found in the section of 
Appendix I indicated: 

Geometric Propertics Seotion 5 
Critioal Pressures Section 3 


Average Critical Pressure Section 3 
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b. Measured Quantities 

1. The length of the colimn was measured with a 6 foot steel tape. 
It is expected that measurements are within .03 inch. 

2. The outer dianeter and thiokmess were measured with a micrometor. 
It is expected that measurements are within „0005 inch. 

Se The radius of zyration was obtainol analytically usin: the outer 
diameter and thiokmess. 

4ο The weight of the column was obtained to the nearest 0.1 ran. 

5. The area was obtained analytically by dividing the tube weicht 
by the product of the length by the density. The tabulated density of 
245-Т3 aluminum alloy is 0.100 pounds per oubio inch. 

6. The tabulated value of the modulus of elastioity of 24S-T3 
aluminium alloy is 10.6 x 105 lbs per inch, This value was verified 
experimentally. 

7. The theoretical critical pressure was computed using a oalcu- 
lator. An excess number of figures was carried to avoid round-off 


errors. 
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